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Abstract 

We study phase segregation in a model alloy undergoing both ordering and decomposition, 
using computer simulations of Kawasaki exchange dynamics on a square lattice. Following 
a quench into the miscibility gap we observe an early stage in which ordering develops while 
the composition remains almost uniform. Then decomposition starts with segregation into 
ordered and disordered phases. The two spherically averaged structure functions, related 
to decomposition and to ordering, were both observed to obey scaling rules in the late 
coarsening stage where the time increase of the characteristic lengths was consistent with 
a(£s -f b). While a was similar for ordering and decomposition at low concentration of 
the minority component, it showed an increase (decrease) with concentration for ordering 
(decomposition). The domain morphology was found to depend on the concentration of 
the minority component, in a way that suggests a wetting of antiphase boundaries in the 
ordered domains by the disordered phase. 
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1. Introduction 

The kinetics of alloy decomposition is a theoretically challenging and technologically im- 
portant problem that has attracted much attention over many years [1-6]. For alloy sys- 
tems where the lattice spacings in matrix and precipitate are the same, and the two-phase 
structure coarsens with time to reduce its interfacial energy, both the morphology and the 
growth rate of the domains seem to have a universal behavior [7]. In such systems, the 
typical domain size grows asymptotically with time t like a power t Q , a = | [8, 9], and the 
essential features of the morphology depend only on the volume fraction / of precipitates 
[9,10]. 

This "universal" morphology is known to break down in cases where the lattice spacings are 
different in matrix and precipitates. The so-induced elastic misfit interactions may change 
the morphology as well as the kinetics of the decomposition. Such systems have been the 
subject of numerous experimental [11-16] and theoretical [17-26] studies in recent years. It 
has been found that anisotropic elastic constants lead to the formation of very anisotropic, 
mostly plate-like precipitates. Nevertheless, their typical size still grows in many cases like 
£3 [24, 15], remaining so even under externally applied stresses [25]. 

An interesting situation occurs in cases where the precipitates are ordered intermetallic 
alloys, such as in Al-Li or in nickel-based superalloys. Here the decomposition (described 
by the conserved order parameter c, that is the local concentration of A atoms in the A-B 
alloy) couples to an ordering transition (described by a non-conserved order parameter). In 
Ni-based superalloys both the ordering of the precipitates and the elastic misfit interactions 
will be present [11-15]. Recent computer work [19, 20, 26, 27] is devoted to the study of 
models that include both effects. While some features of the morphology, like the cuboidal 
shape of the precipitates and their predominant alignment along certain crystallographic 
directions, are clearly due to the elastic misfit interactions, some other features are probably 
due to the ordering of the precipitates. Indeed, it is usually found in these alloys that close 
neighboring precipitates do not merge even though this would reduce the total surface 
energy [11-15]. A possible explanation, proposed in [20], is that these precipitates belong to 
different variants of the same ordered structure so that their merging would lead initially to 
an anti-phase boundary (APB), whose energy is higher than that of the interfaces between 
the ordered precipitates and the disordered matrix. 

In order to address this issue and to gain a better understanding of the influence of atomic 
ordering on the coarsening kinetics and morphology, we have carried out computer sim- 
ulations of an Ising model on a square lattice with nearest and next-nearest neighbor 
interactions, which does not include elastic misfit interactions. There will now be two 
variants of the ordered phase depending on which sublattice is occupied by which type of 
atoms and if a wetting of APB's by the disordered phase reduces the total energy, then the 
domain morphology will be very different from the case where none of the phases (minority 
or majority) is ordered. In particular at larger volume fractions / of the ordered phase, 
e.g. when / > 0.5, APB's are to be expected between the differently ordered domains. It 
is unclear whether the length scale of precipitates in such a system will still grow like £s 
or will behave like t 1 / 2 , the typical growth law for nonconserved quantities. 
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The model we investigate was studied previously by Salmi et al. [28]. Their conclusions 
were mainly qualitative. In the present work we extend their simulations and obtain 
more quantitative information about the joint kinetics of ordering and segregation. We 
concentrate on comparing the growth kinetics and morphology at / < 0.5, / = 0.5 and / > 
0.5 with those obtained for an alloy where none of the phases has an ordered intermetallic 
structure. While the 1 3 law appears to hold in all cases, the coefficient a has very different 
dependence on concentration for the ordering and segregation, the former growing much 
faster for higher concentrations. We also find, in the case of high volume fraction of 
ordered precipitates, clear evidence of the wetting of APB's. This may explain some of the 
results in Al-Li systems [29, 30] where ordered domains with little lattice mismatch are 
formed and supports the interpretation of non-merging precipitates given in [29, 30]. We 
plan to extend our work to include elastic interactions, present in Ni-based alloys [11-15], 
mentioned earlier. 



2. The Model 

We consider an Ising model on a 2D square lattice with nearest neighbor interactions J 
of antiferromagnetic type and next nearest neighbor interactions of ferromagnetic type, 
which in our simulations we took to be — J/2. The total energy U of this system, wrapped 
on a torus containing N = L 2 sites, is 

U = J ~ \ J Yl ai(Tl '' 3 > °' ( 1 ) 

<U'}+ (M')x 

where o\ = 1 if there is an A-atom and a\ = — 1 if there is a B-atom, on the lattice site 
1 = (li, I2) and the first and the second sums are over all nearest and next nearest neighbor 
pairs, respectively. The average concentration c of minority B-atoms is 

c= (I-07/2 

where 

1 

The phase diagram of this system obtained from Monte-Carlo simulation [31, 32] is shown 
schematically in Fig. 1. The critical point for the ordering transition is at c c = 0.5 and 
T c ~ 3.8 J/k and there is a tricritical point at c t ~ 0.27 and T t ~ 1.3 J/k. Above T t 
the transition between the disordered A-rich phase and the ordered phase is of 2nd order, 
while below T the transition is of the 1st order. Fig.l shows the miscibility gap between 
the ordered phase (with practically stoichiometric A-B composition) and the disordered 
phase. 

Our simulations were performed on a square lattice, 128 x 128, with periodic boundary 
conditions, at T = 0.9 J/k at three values for the composition, c = 0.15, 0.25 and 0.35. 
These are all well within the two-phase region of Fig. 1 between the ordered intermetallic 
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alloy with stoichiometric composition AB and the disordered A-rich phase; correspond- 
ing approximately to volume fractions the respective volume fractions of ordered phase 
will therefore be about / = 0.3, 0.5 and 0.7, respectively. The dynamics of the order- 
ing/decomposition process, which are our main concern, were followed using Kawasaki 
dynamics with the Metropolis rule, that is by choosing a nearest neighbor pair at random 
and then exchanging the atoms with probability 

w = min[l, exp(— AU/kT)] 

where AU is the change in the total energy U given by eq.(l). 

In our simulations we start the realizations (computer runs) from randomly gener- 
ated "infinite" temperature configurations. Time is then measured in Monte-Carlo Steps 
(MCS), that is the number of attempted nearest neighbor exchanges per lattice site. 

3. Characterization of Configurations 

Since there are two sublattices within the ordered " antiferromagnetic" structure on 
the square lattice, there will also be two possible variants for the ordered phase, one where 
the A-atoms are on the even sublattice (which we call for short the even variant) and 
one where the A-atoms are on the odd sublattice (odd variant). To represent "snap-shot" 
pictures of the configurations in a way that clearly differentiates between the two variant 
ordered phases as well as the disordered phase, we use the following variable, 

1. = ^i 8 + |E^ ( 2 ) 

where the sum is over nearest neighboring sites, and 

a? = (-l)'i+feoi 

is the staggered spin. The variable 771 takes the value —1 at sites within ordered domains 
of the odd variant and +1 on sites within ordered domains of the even variant. In the 
disordered phase, where almost all atoms are of the same kind, r/i = 0. Intermediate 
values 7/1 = ±|, ±|, ±| appear on the interfaces between ordered and disordered domains 
and at isolated B-atoms dissolved in the A-rich (disordered) phase. Translating this into 
a grey-scale, the two variants of the ordered phases appear black and white, respectively, 
while the disordered phase appears in grey, see Fig. 2. 

The staggered spin local average, 771, can be interpreted as local order parameter for 
ordering and its absolute value averaged over the lattice, 

Vsr = \v\ = \vi\ (3) 

can be considered as a short range order parameter (for the atomic ordering, as opposed to 
decomposition) . It represents the relative amount of ordered phase, irrespective of whether 
the order corresponds to the odd or the even variant. 
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The average over the entire lattice of the staggered spin local average, 

i i 

is equal to the global average of the staggered spin variable and can be considered as a 
long range order parameter. It corresponds to twice the difference in the concentration of 
A-atoms between the odd and the even sublattice. Similar order parameters may also be 
defined for decomposition, although their interpretation is more ambiguous now, 

Psr = ^^2\pi\, p Ir = — J]pi = l-2c (5) 

where 

!/ 1 v-^ \ 
P\ = 2^+4 Z^ ai ')- 

From these definitions it is apparent that the long-range order parameter pi r for decom- 
position is conserved (since c is fixed) while r\i r is non-conserved. One may further notice 
that an alternative way to write r] sr is 

1 1 . 

T] sr = 1 - p sr = - - — ^ = 2 C ( 1 ~ C ) ( 1 ~ «1 ) • ( 6 ) 

(1,1') + 

This expression shows that i] sr is related linearly to the first Warren-Cowley short-range 
order parameter a± generally used in crystallography [33], 

a 1 = l-p AB /c, 

where p AB is the density of B-atoms at nearest neighbor sites of an A-atom. 
To characterize the decomposition kinetics, we use the structure function 

^k=^|^exp(zk.l) ( 7 1 | 2 , (7) 
1 

where k = (fci,/^) is a reciprocal lattice vector, with kj = lixKjjL and Kj = 1,2,..., L 
for j = 1,2. For the system undergoing both ordering and decomposition, the structure 
function (shown in Fig. 3 by means of brightness) is localizing in two areas, around the 
point (it, 7r) (related to ordering) and around the point (0, 0) (related to decomposition) of 
reciprocal space. Note, that ^(o.o) = ^ Pi r an< ^ ^ na ^ = Nr]f r . 

To study kinetics and scaling properties of the two processes we divide the (square shaped) 
Brillouin zone in reciprocal space in two equal areas (bold lines in Fig. 3) : the square for 
ordering, 

O w = {k : |fci - ?r| + \k 2 - ?r| < tt} (8) 
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(that is inside that boundary in Fig. 3) and the remaining four quarters of the square 
for decomposition. Using the periodicity of the reciprocal lattice, the latter can be put 
together as a square, 

O = {k : |fci| + \k 2 \ < tt} (9) 

shown separately in Fig. 3 with higher resolution for values of the structure function. To 
test scaling behavior of the structure function, we sphericalize it within the two squares 
fij, related to ordering (y = it) and to decomposition (y = 0), separately. We define: 

s v (k) = s *i E ^ ( 10 ) 

A v (k) A„(k) 

where 

A v (k) = {k' : k - 5k < |k' - (u, v)\ < k + 5k} 

is a ring of width 25k and median radius k = |k| = {k\ + k%)^ , centered at (u, v). The 
sphericalized structure function was actually computed for the values k = 2irK/L, K being 
an integer: K = 1,2, [L/2y/2\, and choosing 5k = n/L. Since the configurations on the 
lattice are evolving with time, we are also using the notation S l/ (k,t) for the sphericalized 
structure function computed for the configuration at time t. 

The scaling hypothesis for the late stages of coarsening, in systems undergoing either 
phase segregation or ordering, is that the system is characterized by a single scale X(t). This 
means in particular that for values of k~ l which are large compared to the lattice spacing 
and small compared to the size of the system S(k,t) should (up to a time dependent factor) 
be a function of only one variable, X(t)k. In our case we need a priory two characteristic 
lengths, \ v (t), v = 0, n. We therefore hypothesize that for the system under consideration 

S v {k,t)fiiA v {t)F v (kX v (t)) as t — > oo, (11) 

where S„(k, t) is the macroscopic, formally the infinite volume limit, for fixed k ^ 0, of the 
corresponding sphericalized structure function at time t [9]: 

S u (k,i) — > S u (k,t) as L — > oo 

We shall show, in section 4, that this generalized scaling hypothesis is indeed consistent 
with our simulation data. 

One of the common ways of defining a characteristic length is to use inverse moments 
of the structure function. Since we are restricted in simulations to systems which are very 
small on the macroscopic scale, great care has to be taken in computing from the available 
data quantities relevant for macroscopic systems. This is doubly so for our system where 
we have both ordering and segregation. We therefore had to find a way which minimizes 
distortions of the time dependence of the \ v due to the finite size of the system. To achieve 
this we defined the moments 

k b k b 

k v (t) =J2k 1+r "S v (k,t)/J2k rv Mk,t) (12) 



6 



with ro = 1 and r n = 2, where the limits of summation, k a and k(,, were themselves 
scaled in accordance with k v (t), by putting k a ^(t) = (1 ± h)k u (t); so at different times the 
moments are calculated from similar intervals. Taking h = .7, the interval around k v was 
large enough to give sufficient statistics for the summation, but still small enough to avoid 
any contribution from the structure function at very small or very large k, which might be 
affected by finite size effects. The k v (t) so obtained appeared to follow the position of the 
maximum of k Tv S v (k,t) in eq. (12). 

With the procedure described above, finite size effects on k u (t) could be minimized. 
Indeed, additional runs performed on a smaller lattice (64x64) gave consistent results 
for the time evolution of k u (t). Of course this still leaves considerable uncertainties in 
determining k u (t). For the size of system we are dealing with, the gaps in the fc's, which 
are equal to 2tt/L, limit any finer resolution independent of the number or duration of the 
runs. 

4. Results 

Snapshot pictures of typical configurations obtained at T = 0.9 J / k for alloy compositions 
c = 0.15, 0.25 and 0.35 with lattice size L = 128 are shown in Fig. 2. The value for the 
locally averaged spin 771, defined in eq.(2) is shown at 125, 8000 and 64000 MCS. Three 
different types of domains appear in these pictures in white, black and grey corresponding 
to ordered phases on the even and the odd sublattice and the disordered A-rich phase, 
respectively. In the case c = 0.25, the grey phase covers about half of the specimen volume, 
while at c = 0.35 it covers less and at c = 0.15 more. This is, as already mentioned, in 
qualitative agreement with the equilibrium phase diagram [28, 31, 32] where the predicted 
volume fractions of black or white (that is ordered) phase are / = 0.3, 0.5 and 0.7 for c = 

0. 15, 0.25 and 0.35, respectively, see Fig. 1. In the case of decomposition without ordering, 
when the miscibility gap is symmetric around the critical composition, the behavior of an 
alloy with / = 0.3 would be identical to the one at / = 0.7, with the two phases exchanged 
in their respective role. In the present case, however, it becomes clear from looking at 
Fig. 2 that there is a completely different behavior on different sides of the tricritical 
composition. At c = .15 the pictures show isolated ordered droplets within a disordered 
matrix, while for c = .35 they show the disordered phase surrounding the ordered domains, 

1. e. we see an ordered alloy with the antiphase boundaries between the two ordered variants 
wetted by the disordered phase. 

4.1. Early Stage Behavior 

We observe that the short range order parameter, r] sr comes rapidly to the steady value 
2c, corresponding to complete ordering. We find that the relative deviation becomes small, 
\i] sr /2c— 1| < 10 -2 , by the time t ~ 10 2 , indicating that local ordering is well developed by 
that time. At the same time a domain structure starts to be visible (see Fig. 2). Ordered 
domains of both kinds (even and odd variant) occur in about equal fractions. Domains of 
pure minority component are practically not observed (in the original spin configuration) 
after time t ~ 10 2 . Then segregation of ordered and majority phases develops. 
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We further observe that the sphericalized structure function for ordering, S n (k), de- 
velops a maximum localized at k ~ by the time t ~ 10 (starting from times t ~ 1) 
which grows sharper as t increases (Fig. 4.1), while the sphericalized structure function 
for decomposition, So(k), develops a maximum at k ~ ko later, starting at t ~ 10 2 (Fig. 
4.2). This shows, in agreement with what was found in [29, 30, 34], that, initially, ordering 
is much faster than decomposition. 

Moreover, the sum of the structure function components related to ordering 

N S7r = J2 S * ( 13 ) 

approaches the constant value 2iVc, which corresponds to the area of ordered phase at 
complete phase separation. Similarly, the sum of the structure function components related 
to decomposition (that is, within the square f^o with exclusion of the origin), 

Ns = Ys S * ( 14 ) 

% 

also reaches a steady value Nsq, with Sq ~ (1 — 2c) 2c, which was to be expected, since 
there is a general relation between sq and s n : 

S7r + s = l- S {0fi) /N = 4c(l-c). (15) 

The relative deviation of sq from its equilibrium value, |so/sq — 1|, becomes smaller than 
lO- 2 , for times t > 10 2 . 

4.2. Scaling Analysis 

As already mentioned, both spherical averages of the structure functions, S n (k), and 
So (k) defined in (10), have a well defined single maximum by the time t ~ 10 3 or 10 4 
(see Figs. 4. a and 4.b). Using these functions to compute the moments defined in (12) we 
find that the two characteristic dimensions related to ordering and decomposition, X u (t), 
satisfy [35, 9] 

\„(t) = 2Tr/k u (t)f*ia u (t*+b u ) (16) 

for times greater than 10 3 . This can be seen in Figs. 5. a and 5.b, for ordering and 
decomposition, respectively. The constants a v and b v are given in Table 1. Equation (16) 
shows that the growth of all characteristic dimensions is consistent with the fimiliar t 1 / 3 
growth law [8, 35] although an unambiguous determination of the growth law exponent 
would only be possible with much longer runs on a considerably larger lattice. 

We now investigate the scaling hypothesis for both functions which requires that 
S v (k, t) can be written in the form 

S v (k,t) » B v k' v {t)~ 2 F v (k/k' v (t)), (17) 
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where we have used the constancy of the sums in (13) and (14); k' v (t) may differ from the 
previously defined k v (t) by a time-independent (but composition-dependent) factor, i.e. 

2n/kl(t) =a' v {t* + b v ). (18) 

A standard way to choose the constants a' and B is to require that the scaling function 
F satisfy the normalization condition 

max[F (x)] = F (l) = 1. (19) 

This permits comparison of characteristic parameters at different concentrations. 

To better locate the maximum of the function Fq(x), we tried various empirical fits 
starting with a one parameter fitting expression proposed for decomposition in two dimen- 
sions, Eq. (14) in [9]. While this works well for pure segregation and gave a reasonable fit 
to our data, it appeared to miss the stronger asymmetry, about the maximum, present in 
our system. This was particularly apparent at c = .35 when the morphology is indeed very 
different from the pure segregation case. We found by trial and error that the following 
formula 

F (x) = ax 4 /(^ + x 4 )[p+(x 2 -l + 8) 2 ] r , a,/?, 7 >0, (21) 

r = 3/4, gave a good fit. The form (21) with 3/4 replaced by (d + l)/4 is a generalization, 
satisfying Porod's law (Fq(x) ~ x~( d+1 ^&t x — > oo) [36], of a heuristic formula suggested 
in [10] for 3D alloys (where d = 3). Definition of the standard form function Fq(x) yields 
two relations between parameters 

a = (l + 7 )(/3 + 5 2 ) r , 7 = rS/((3 - rS + 5 2 ) 

so, two of the four parameters, preferably (3 and 5, are independent. These were fitted 
from the data. 

For F-X the situation was more difficult, because there was no scaling close to (it,tt). 
The reason for this is that S(tt, tt) is related to the (non-conserved) long-range order pa- 
rameter (S(tt, tt) = Nr)f r ), which depends strongly on how the ordered phase is distributed 
over the two sublattices. Further away from the (tt, 7r)-point, the structure function de- 
scribes the short-range order, which is not sensitive to the difference between the two 
sublattices. Because of that, a' n and were determined by fitting of F 7r (x) in the region 
where scaling holds, to the square of the Lorentzian function, 

F n (x) = l/(l + x 2 ) 2 , (22) 

As shown in Fig. 6, we observe scaling of both structure functions (for ordering and for 
decomposition) for times t ~ 10 through t ~ 10 5 . The obtained values for the constants 
P, 5, a! v and B v are given in Table 2. We find that a'^, giving the rate of ordering kinetics, 
increases while the rate a' of decomposition (coarsening) decreases with the amount of 
ordered phase at equilibrium. This is presumably due to the very different morphologies 
at high and low concentrations discussed earlier. 
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Finally, the dependence of full-width at half maximum (FWHM) of the scaling func- 
tion on the volume fraction of B-rich phase is given in Fig. 7 where it is compared with 
the FWHMf ei .ro obtained from simulations on the two-dimensional Ising model with at- 
tractive nearest neighbor interaction between like atoms [9] . The FWHMf erro is practically 
temperature independent (for 0.34 < T/T c < 0.8, T c being the critical temperature) and 
symmetric around the volume fraction / = 0.5. In the present model, however, the scaling 
function is wider in all cases, indicating a smaller amount of regularity in the positions of 
the precipitates. Moreover, there is a clear asymmetry between / as 0.3 and / ~ 0.7, which 
is not surprising since the morphologies of the two-phase mixtures are also quite different 
at these compositions (cf. Fig. 2). 



Discussion 

The general picture emerging for the ordering/decomposition process in the present model 
is that the phase transformation starts with the development of (short-range) order followed 
by the appearance of well-ordered domains inside a disordered matrix. At this stage, the 
ordering (as characterized by the short-range order parameter ?y sr , eq.3) has reached its 
maximum value and the further development corresponds to the coarsening of domains 
by long-range diffusion, much like in a situation where decomposition into two pure (or 
disordered) phases occurs. The fact that ordering appears prior to decomposition has been 
observed experimentally [29, 30, 11-15] and can be explained by the fact that the diffusion 
of B-atoms is required only over short distances for the development of atomic order but 
needs to cover long distances for the decomposition into phases with different amounts of 
B-atoms. 

The coarsening kinetics of the domains, measured by the correlation lengths for decompo- 
sition, , and for ordering, fc" 1 , is consistent with a growth of domain size proportional 
to t 1 / 3 , in agreement with experiments on ordered precipitates [29, 39, 11-15]. This growth 
law, as well as the observed time-scaling of the structure function, are a signature of the 
coarsening of two-phase structures under the influence of interfacial energy [9, 35, 37]. The 
most striking difference between the present case and a conventional coarsening process 
(with two disordered or pure phases) is the composition dependence of the morphology. 

In particular, the lack in symmetry between the case c = 0.15 (with volume fraction 
/ ~ 0.3) and c = 0.35 (with volume fraction / m 0.7) may be related to the presence of 
ordered domains of two kinds (on the two sublattices) with some additional positive energy 
appearing (as pointed out, e.g., in [29, 30, 34]) on their interface when they are coming in 
contact due to their growth (anti-phase boundary). These anti-phase boundaries (APB's) 
appear to be wetted by the disordered phase which reduces the overall energy of the sys- 
tem. Hence, even though the ordered phase is the majority at c = 0.35, individual ordered 
droplets do not join up but stay separated by narrow channels of disordered phase (see 
Fig. 2c). The composition-dependent morphology may also be at the origin of the dramat- 
ically different behavior of the typical length scales for ordering and for decomposition 
Aq, when the alloy composition is changed. 
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Table 1. Parameters a v and b v characterizing the growth (16) of of the characteristic 
lengths defined in (12), at different concentrations c of minority component. 



c a n bn a b 



0.15 


0.78 


9.0 


0.77 


7.0 


0.25 


0.94 


11.0 


0.67 


5.5 


0.35 


1.7 


10.0 


0.475 


9.0 



Table 2. Parameters a' v and B u used for the scaling of the structure functions, eqs.(17) 
and (18), and parameters j3 and S used for the fitting of F (x) with eq.(21), as a function 
of the concentration of minority atoms, c. 



c 


< 




a' 


Bo 


P 


5 


.15 


0.82 


3.1 


0.97 


0.53 


1.47 


0.13 


.25 


1.01 


5.1 


0.79 


0.75 


0.98 


0.065 


.35 


1.75 


6.7 


0.65 


0.37 


3.5 


0.55 
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FIGURE CAPTIONS 



Fig. 1. Schematic equilibrium phase diagram (concentration of B-atoms, c, versus temper- 
ature, T) for the model A-B alloy. The circles show the conditions at which our computer 
experiments were performed. 

Fig. 2. Plots of configurations of the local order parameter, 771, (by means of brightness 
i]i = —1 being white and rji = +1 being black), at three time points, frs = 5, 20, 40 (time, 
t in units of MCS, from top to bottom) for each of three concentrations, c = 0.15 , 0.25, 
and 0.35. The disordered phase has r]\ and appears grey. The lattice size is 128 x 128. 

Fig. 3. Plot (by means of brightness, black meaning high values) of the structure function, 
Sk, in the entire Brillouin zone (top) and in its parts related to decomposition (bottom, 
with higher contrast) at = 5 (time, t is given in units of MCS), for the concentration 
c = 0.25. 

Fig. 4. (a) Plots of the spherically averaged structure function for ordering, S 7T (k), at 
= 2, 5, 10 (time, t is given in units of MCS) for concentration c = 0.25. (b) So(k) at 
times t 1 / 3 = 5, 10, 20. 

Fig. 5. Dependence of the characteristic wavelength for (a) ordering, A^, and (b) decom- 
position, Ao, on time, t (in units of MCS) at three concentrations: c = 0.15 (diamonds), 
0.25 (stars) and 0.35 (triangles). 

Fig. 6 (a) Scaling of the structure function spherical average for ordering, F n (k/k 7T ) where 
k n = a'^(t l i + b n ), for c = 0.15,0.25 and 0.35. The full line is a fit with eq.(22). (b) 
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Scaling of the structure function spherical average for decomposition, F (k/k o ), where 
ko = «o(t3 + bo) (at the same three concentrations). The full line is a fit with eq. (21), 
with the parameters given in Table 2. 

Fig. 7. Full-width at half maximum of the scaling function for decomposition normalized 
according to eq.(19), as a function of the volume fraction / of the B-rich phase (crosses). 
Data obtained with a two-dimensional Ising model with attractive nearest neighbor inter- 
action between like atoms (taken from [9]) are also shown for comparison (diamonds). 
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